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We provide a simple method for the calculation of the terms c„ in the Zassenhaus product e""'"'' = 
e" • e*" ■ ri^2 ^'^^ non-commuting a and b. This method has been implemented in a computer program. Fur- 
thermore, we formulate a conjecture on how to translate these results into nested commutators. This conjecture 
was checked up to order n = 17 using a computer. 
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CN . I. INTRODUCTION 

The product of the exponentials of two non-commuting variables x and y may be expressed in terms of the Baker-Campbell- 
Hausdorff (BCH) series 

^ ; e^e'' = e^+''+^"=2^". (1) 

T— I . The terms Zn of the sum may be written as linear combinations of words W of length n consisting of letters x and y, 
> . 

\D] = J2 <W)W{s^,...,Sr,), (2) 

^— I , Wisi s„) 

O ■ 

^f) . where each word W{si, . . . , s„) is a product of n factors = a; or s; = y. The sum in Eq. (|3 is over all possible different 
words, i.e. the sum over W has in principle 2" terms. Some of these terms vanish, because the corresponding coefficient z{W) 

, equals zero. The coefficients z{W) may be determined in various ways. A simple method to determine these coefficients has 

' been suggested recently by Reinsch Q|. His method is easily implemented in a computer program yj^jS]. Another method was 

(~| presented by Goldberg 0], and his method was implemented in a program by Newman and Thompson 

0^' The terms z„ of the BCH series may be expressed as linear combinations of nested commutators of x and y. This was originally 

Jh , shown by Baker, Campbell, and Hausdorff Q]. However, explicit determination of the coefficients was difficult for general n. 

. Dynkin @,0] showed that the coefficients are in fact easily obtained, so that the z„ can be represented as 

a ■ 

6: ^» = - E z{w)[[...[[s,,s2],s3],...],s„] (3) 

> 



X 



n 

W{si, 



where [[. . . [[si, S2], S3], . . .], s„] is a direct translation of the word W into a left normal nested commutator, i.e. the order of 
the letters in the commutator and the word are the same. The representation of the z„ in terms of commutators is not unique 
' due to the Jacobi identity [[a:,?/j,z] + [[y, z],^] + = and similar identities for higher order commutators. Some 

of them are discussed by Oteo psi]. Oteo also formulates a conjecture concerning another translation of the BCH terms into a 
linear combination of commutators. This translation consists of fewer terms than Dynkin's translation, but a general proof of its 
validity is not known to us. We have checked the validity of Oteo's conjecture up to order n = 17 on a computer 
It was first shown by Zassenhaus |0] that there exists a formula analogous to the BCH formula for the exponential of the sum of 
two non-commuting variables a and b, 



ga+6 ^ ga . gb 



He^", (4) 



n=2 



which is known as the Zassenhaus product formula. The individual terms c„ will be called the Zassenhaus exponents in the 
following. They may also be written in terms of words of length n consisting of the letters a and b, i.e. 

Cn = E ciW)Witi,...,tn) (5) 
W(ti,...,t„) 
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with ti = a or ti = b. It is the goal of this note to provide a simple method for the calculation of the coefficients c{W) as well 
as to propose a suitable computer implementation. Our method has been developed in analogy to the procedure proposed by 
Reinsch (H for the BCH terms. 

The Zassenhaus exponents may also be obtained in terms of nested commutators as shown e.g. in Refs. Dynkin's 
theorem 0] provides a (generally valid) translation of words into nested commutators, if such a translation exists. Therefore, 
it is possible to directly translate the expressions for the Zassenhaus exponents (|5} into a linear combination of left normal 
nested commutators as in Eq. (|3}- In analogy to the conjecture by Oteo for the BCH terms, we formulate here a conjecture 
concerning another translation of words into left normal nested commutators for the Zassenhaus exponents involving fewer 
terms than Dynkin's translation. Using a computer, we found this conjecture to be valid for the Zassenhaus exponents up to 
order n = 17, but at this time we cannot provide a general proof. 



11. THE ZASSENHAUS EXPONENTS 



In this section we state a corollary which allows a recursive determination of the Zassenhaus exponents. 

Let Ti, . . . , r„ be commuting variables. In terms of these variables we define three upper triangular (?? + 1) x {ji + 1) matrices 
H, K, and L with matrix elements given by 



k 



i-iy_ 

k—i 



for l<i<j<n+l and zero otherwise. These matrices may be expressed as exponentials of the (n + 1) x (n + 1) matrices 
P and Q defined by 

Pij — Qij = ^i+l,jTi (8) 

where Sij is the Kronecker symbol, 

H = cxp{P + Q), K = cxp{-P), L = cxp(-Q). (9) 
Furthermore, we define an operator U which operates on products p of the variables r, 

with G {O,l}fori ~ l,..,n. The operator L/ "translates" such a product p into a word consisting of letters a and b according 
to the following rule: If = 0, rf' is replaced by an a, and if fii = 1, rf ' is replaced by a b. The index i determines the 
position of the letter in the word. E.g., for n = 6 the product p = '''i ''"2 '''3'^4'''5 '''e — '''i'''3'^4 is translated as follows: 

U{p) ^ U{tIt^t^tIt^t^) = [/(riT3T4) = babbaa. 

The operator t/ is a vector-space isomorphism mapping the space of polynomials in the r-variables (with /i; = or /i^ = 1) into 
the space of words of length n. 

Corollary 1 The Zassenhaus exponent C2 defined in Eq. is obtained in terms of the 3x3 matrices L, K, H as C2 = 
U{L ■ K ■ iJ)i,3. For n > 2, the Zassenhaus exponents c„ are given in terms of the corresponding (n + 1) x (n+ 1) matrices as 



t/((e 



-'-'■■■■■^-"'■L-K-H)^^ ) (10) 



Here, Cm (1 < m < n) are the Zassenhaus exponents written in terms of the (71 + 1) x (n + 1) matrices P and Q, and the index 
1, 71 + 1 indicates the upper right element of a matrix. 

This corollary permits a recursive determination of the Zassenhaus exponents. In fact, due to the special structure of the matrices 
P and Q all exponentials in Eq. ( I10> are obtained as finite sums, and the whole calculation can be done in a finite amount of 
steps either by hand or on a computer A suitable computer implementation will be presented in section lTvl 



A. Examples 



Before proving Corollary 1 we work out a number of examples. 
For n = 2 we need to use the 3x3 matrices given by 



1 -Tl ^TlT2 
L = I 1 -T2 \ , K 

1 

1 (1+Tl) i(l+Tl)(l+r2] 

H = \ 1 (l + rs) 

1 



1 -1 I 
1 -1 

1 



Then, 



LKH = 



1 in - ir2 

1 




and the second Zassenhaus exponent takes the form 



C2 = U{L-K-H)i 



n+l 



{ha — ab). 



For n = 3 we need to use the 4x4 matrices 



L 



K 



H 



/I 


-Tl 


1-^ ^ 1 
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-T2 








1 










/I 


-1 


2 6 





1 










1 -1 







1 1 
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(1 + 


ri) + 
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T1T2T3 
■T2T3 

1 



(1 + ^2) 

1 





1=1 

1=2 

1 



and the matrices P and Q defined in Eq. (|8} 



10 

10 

1 





and Q = 



















T2 














T3 


u 












It follows from Eq. (tTll that 



-C2 



exp 

(l i(T2-Ti) 

1 

1 

\0 





i(r3-r2) 

1 
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Therefore, the third Zassenhaus exponent takes the form 



C3 = C/(e-^^.L./i.ff)^„^^ 

^ /2 11 1 11 

= U (^-nra - -r, - -r,r, -r,r, + -r, + -^r, 

2 11111 

= —hab aba abb bba H — baa H — aab. 

3 3 3 3 6 6 



In an analogous way one obtains for C4 



C4 = ——aaab + —aaba + —aabb — —abaa 
24 8 8 8 

abab abbb H baaa H — baba 

4 8 24 4 

3 13 1 

+—babb bbaa bbab H — bbba. 

8 8 8 8 



These results all agree with the standard results given in the literature, see e.g. Ref. fllll . 



III. PROOF OF COROLLARY 1 

The (71 + 1) X (?i + 1) matrices P and Q have non-zero elements only in their first superdiagonal. A product of m factors P or 
Q is a (n + 1) x (rt + 1) matrix which contains non-zero elements only in the jnth superdiagonal. In particular, a product of 
n factors P or Q has only one non-zero element in its upper right corner A product of k factors P or Q with fc > 71 is a null 
matrix. 

Each Zassenhaus exponent c„ is a linear combination of words of length n. As a consequence, C„ is a (n + 1) x (n + 1) matrix, 
which has a non-zero entry only in its upper right corner, and for each n e N 



e^+« = e^ • e« 



Yle"^ (12) 



is obtained in terms of a finite product with n factors. Therefore, one obtains 

e'^" - e-C"- ..... e-c^ . e-^ . e"^ . e^+'3. (13) 
The exponentials can be calculated from a finite sum. 

Since C„ only has one non-zero element in its upper right corner, it holds that e*"" — I + Cn, and we find 

= (e-'^"-.....e-'^^.e-Q.e-^.e^+Q)^„^^, (14) 



where Ai = P ot Ai = Q, and the sum runs over all different matrix products Ai . . . An- The next step is to show that the 
upper-right element of a product matrix Ai . . . An is given by a product of r^, and that the indices on the variables determine 
the positions of the Q's in the matrix product Ai . . . An- This has been shown in Ref. 01 ^nd is not repeated here. Applying the 
operator U on this product of then transforms the result Eq. ( I14> into a linear combination of words in terms of the letters a 
and b- This proves Corollary 1 . 



IV. COMPUTER IMPLEMENTATION 



The following Mathematica program implements Corollary 1 . Calling it will return the Zassenhaus exponent c„ in terms of the 
variables a and b- 

The program consists of three parts: First the matrices L, H, P and Q are defined. Then the product of exponentials as 
required by Corollary 1 is calculated (starting from n — 2), and finally the translation U is implemented. The program works 
with strings in order to prevent Mathematica from sorting the words alphabetically. 
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ZH[n_,a_,b_] := Module [{ C, L, K, H, P , Q, m, t , r, i , j , k, u, z } , 

C[2] = (t[l]-2 t [2] ) /2- (t [1] t[2]"2)/2; 

For [m =2, m <= n, m++, 

L = Table [ (-1) ' (i+j) / (j-i) ! Product [t [ k] ,{ k, i , j-1 }],{ i , m+1 },{ j , m+1 }] ; 
K = Table [ (-1) ' (i+j) / (j-i) ! , { i , m+1 } , { j , m+1 } ] ; 

H = Table [1/ ( j-i) ! Product [( 1+t [ k] ),{ k, i , j-1 }],{ i , m+1 },{ j , m+1 }] ; 

P = Table [KroneckerDelta [i + 1, j ],{ i, m+1 },{ j , m+1 }] ; 

Q = Table [KroneckerDelta [i + 1, j] t [ i ] , { i , m+1 } , { j , m+1 } ] ; 

C [m] = Expand [ 

((Dot @(a Table [MatrixExp [-Sum [ 

r = (List @(a C[m-u][[z]]) /. {t[i_] -> P, t[i_]"2 -> Q}; 
r[[l]] (Dot (?@ Take [r, -Length [r] +1] ) , 
{ z, Length [C [m-u] ] } ] ] , 

{u, l,m-2}] ) .L.K.H) [ [1, m+1]] Product [t [ j ],{ j , m} ]]] ; 

Sum[r = (List (5(3 C [n] [ [k] ] ) 

/■ {t[i_] -> ToString[a], t[i_]"2 -> ToString [b] } ; 
r[[l]] (StringJoin @Q Take[r,-n]), {k, Length [C [n] ] } ] ]; 

More elegant but less readable Mathematica implementations than the one given above are possible, e.g. using Nest List 
instead of a For loop. Since the Zassenhaus exponents for larger n are rather lengthy expressions, the computer needs a 
significant amount of memory for this calculation. On a standard personal computer with 2 GB of memory we could obtain the 
Zassenhaus exponents up to n = 17 within about one hour of computer time. 



V. EXPRESSION IN TERMS OF COIVIMUTATORS 



In the introduction we briefly discussed Dynkin's translation la, LL of words into commutators, which is applicable, whenever 
such a translation exists. Since it is known that a representation in terms of commutators exists for the Zassenhaus exponents 0, 
[Toll , we may directly use Dynkin's prescription in order to obtain an explicit representation of the Zassenhaus exponents in terms 
of commutators. We checked the validity of this procedure using a computer up to order rt = 17. A translation for the BCH 
terms into commutators involving fewer terms than Dynkin's prescription was proposed by Oteo 0]. To our knowledge the 
validity of this translation has never been proved in general. Oteo showed it to be valid up to order n = 10, and we checked this 
conjecture up to order n = 17 on a computer 

In analogy to Oteo's prescription for the BCH terms, we now write down an expression for the Zassenhaus exponents 

c„ = c{W)W{tu..,tn) (15) 

w(ti,..,t„) 

in terms of commutators. The words W consist of letters a and b; na{W) counts the number of a's in that word. Analogously, 
nb{W) = n — na{W) counts the number of 6's in the word W . We conjecture that the Zassenhaus exponent c„ may be expressed 
in terms of the left normal commutator as follows 

^" = E (16) 

»-(tl,-.,t„) ' 

Here we only sum over words starting with the letters ha. Similarly, one may write 

c{W) 



E 



na(W) 

t, =a, t, = b 



[[tiMMi-]M (17) 



and only sum over words starting with the letters ah. 

We checked this conjecture up to order n = 17 using Mathematica and compared results up to order n = 6 with expressions 
given in the literature, e.g. in Ref. II lH . 
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A. Example 

For n — 4 one finds the following representation of the Zassenhaus exponent in terms of words 

C4 = ——aaab + —aaba + —aabb — —abaa 
24 8 8 8 

——abab — —abbb + —baaa + —baba 
4 8 24 4 

3 13 1 

~\-—babb bbaa bbab H — bbba. 

8 8 8 8 

According to our conjecture this may be translated into nested commutators as 

C4 = ^ma],a],a] + ^[[[b,a],b],a] + ^[[[b,a],b],b] 
= -^[[[a,b],a],a] - ^[[[a,b],a],b] - i[[[a, 6], 5], 5]. 
This result agrees with results given in the hterature (e.g. in Ref. [10]). 



VI. CONCLUSION 

We developed a method for the calculation of the Zassenhaus exponents c„ in the Zassenhaus product formula e'^+^ = e° ■ e** ■ 
J|J^2 s^" for non-commuting a and b. The method is given in Corollary 1, from which we obtain the Zassenhaus exponents 
in terms of words. It appears that our method is simpler and faster than previous methods (see e.g. Ref. lHoll ). We provide a 
suitable computer implementation. 

We would like to mention that the method presented here can be easily generalized to the Zassenhaus product formula for q- 
deformed exponentials. (The Zassenhaus formula for g-deformed exponentials is discussed e.g. in Ref. fllll and references 
therein.) 

Furthermore, we formulated a conjecture on how to translate the Zassenhaus exponents given in terms of words into a form 
in terms of left normal nested commutators. This representation involves fewer terms than a translation based on Dynkin's 
theorem Q] and has been found to be valid for Zassenhaus exponents up to order n — 17 using a computer We expect a 
proof of our conjecture to be possible along the lines of Dynkin's proof for his representation of the BCH formula in terms of 
commutators. This proof shows essentially by direct calculation that the conjectured commutator representation is equivalent to 
the representation in terms of words. We will address this issue in a forthcoming paper. 
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